# read in spatial data
counties <- read_sf(here("_posts", "2021-03-14-map", "ca_counties", layer = "CA_Counties_TIGER2016.shp")) %>%
clean_names()
oil_spills <- read_sf(here("_posts", "2021-03-14-map","Oil_Spill_Incidents", layer = "Oil_Spill_Incident_Tracking_%5Bds394%5D.shp")) %>%
clean_names()
# check crs
# st_crs(counties) # WGS 84
# st_crs(oil_spills) # WGS 84
# set to interactive viewing mode
tmap_mode("view")
# make exploratory map
tm_shape(counties)+
tm_borders("black")+
tm_shape(oil_spills)+
tm_dots("red")